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The origin of the suppression of the longitudinal spin Seebeck effect by applied magnetic fields 
is studied. We perform numerical simulations of the stochastic Landau-Lifshitz-Gilbert equation of 
motion for an atomistic spin model and calculate the magnon accumulation in linear temperature 
gradients for different strengths of applied magnetic fields and different length scales of the temper¬ 
ature gradient. We observe a decrease of the magnon accumulation with increasing magnetic field 
and we reveal that the origin of this effect is a field dependent change of the frequency distribution 
of the propagating magnons. With increasing field the magnonic spin currents are reduced due to a 
suppression of parts of the frequency spectrum. By comparison with measurements of the magnetic 
field dependent longitudinal spin Seebeck effect in YIG thin films with various thicknesses, we find 
that our model describes the experimental data very well, demonstrating the importance of this 
effect for experimental systems. 

PACS numbers: 75.30.Ds, 75.40.Mg, 75.76.+j 


Spin caloritronics is an emerging research field promis¬ 
ing spintronic devices with new functionalite^ which rest 
on the combined transport of heat and spin [l|, ■ A key 

tool that stimulated this field is the inverse spin Hall ef¬ 
fect 1^ , which allows one to detect pure spin currents 
due to a transformation into measurable charge currents 
via spin orbit coupling. Using this indirect measurement 
technique Uchida et al. measured the so-called spin See¬ 
beck effect (SSE) [ 1 ]. The measurements showed that in 
the ferromagnetic insulator YIG a pure spin current is 
created due to an applied temperature gradient. 

These findings t rigg ered a variety of further studies in 
different groups juMlSf to investigate the origin of the de¬ 
tected spin currents. Interface effects due to a proximity 
effect which creates a magnetization in the normal metal 
were discussed [l^ . but various groups showed that the 
effect in the YIG/Pt wstem cannot be responsible for 
the measured signals 0 i, E A first theoretical 
description of the magnonic spin Seebeck effect was de¬ 
veloped by Xiao et al. [I3- With a two temperature 
model including the local magnon and phonon temper¬ 
atures the measured spin Seebeck voltage was shown to 
originate from the local difference between magnon and 
phonon temperature at the Pt interface caused by ther¬ 
mal spin pumping. Later it was shown by spin model 
simulations that such a temperature difference can be 
caused by magnons traveling from the hotter towards the 
colder part of the sample [l8| . By studying the thickness 
dependence of the SSE for thin YIG films, Kehlberger et 
al. Q found an increasing spin Seebeck signal saturat¬ 
ing above a critical length scale. This length scale can be 
referred to as the mean magnon propagation length. Sim¬ 
ilarly Agrawal et al. measured the time evolution of the 
SSE in the sub-microsecond regime and explained their 
time dependent increase of the SSE by a dependence on 
the characteristic length scale of magnon propagation of 
a few hundred nm 121. 


Further theoretical studies focused on the magnonic 
origin of the SSE [l^, also predicting thickness de¬ 
pendent effects. Hoffman et al. [19j studied the SSE an¬ 
alytically within the framework of the stochastic Landau- 
Lifshitz-Gilbert equation. They assume that mainly 
magnons with hu k, k^T contribute to the effect, but 
Boona et al. proposed that also sub-thermal magnons 
with lower energy can contribute or even dominate 211. 
To reveal, which magnons contribute to the effect, field 
- dependent measurements that tailor the magnon spec¬ 
trum are a key. First such measurements were made 
available in a recent publication by Kikkawa et al. Q 
who studied the reduction of the SSE in the regime of 
high magnetic fields and report a suppression of the mea¬ 
sured SSE signal. However no theoretical explanation 
was given that reproduces the results completely and al¬ 
lows for a direct comparison of experiment and theory. 

In this paper, we study the influence of external mag¬ 
netic fields on the SSE. Applying external fields ma¬ 
nipulates the frequency distribution of the propagating 
magnons by increasing the frequency gap Wmin of the 
frequency spectrum. Therefore, the number of propa¬ 
gating magnons can be reduced by increasing the exter¬ 
nal magnetic field. We present numerical simulations of 
the magnon accumulation in linear temperature gradi¬ 
ents and its dependence on external magnetic fields. For 
this study an atomistic spin model is used, which provides 
realistic spin wave spectra beyond parabolic approxima¬ 
tion and without any artificial cut-off due to discretiza¬ 
tion effects. To understand the underlying physics, an 
analytical model is derived that explains the effect of the 
magnetic field within the framework of linear spin wave 
theory. Finally, our theoretical work is directly com¬ 
pared with field dependent SSE measurements of YIG 
thin films with different thicknesses, which cover a large 
range across the critical length scale showing good agree¬ 
ment. 
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For the numerical simulations we use an atomistic spin 
system consisting of localized spins with normalized mag¬ 
netic moment on a cubic lattice. Our model 

Hamiltonian 'H includes Heisenberg exchange interaction 
for nearest neighbors with exchange constant J, a uni¬ 
axial anisotropy with an easy-axis in z-direction with 
anisotropy constant dz, and an external field B = 
parallel to the easy axis, 

n = -^J2 E ■ (1) 

<i,3> i i 

The dynamics of each single spin are described by the 
stochastic Landau-Lifshitz-Gilbert (LLG) equation, 

consisting of a precession around its effective field and 
a phenomenological damping with damping constant a 
[22L [ 23 j. 7 is the gyromagnetic ratio and the effective 
field Hi is given by 

f)qj 

H.=-^ + C,W- ( 3 ) 

The temperature is included as additional noise term 
Ci(t) in the effective field Hi fulfilling 

m) = 0 (C(o)cj(t)) = (4) 

The dynamics of the systems are studied by numerical in¬ 
tegration of these equations using the Heun-method [2^ . 

In the systems considered, the temperature is spatially 
dependent, including a linear temperature gradient in z- 
direction with a constant slope over a distance L and 
two heat baths at the two ends of the system as shown in 
FIG. [TJ The dimension of both heat baths is chosen to be 
large enough to minimize finite size effects in the area of 
the gradient. The given temperature profile describes the 
temperature of the phononic heat bath and it is assumed 
that the phonon temperature remains constant during 
the simulation. 

We study the manipulation of the magnon accumula¬ 
tion by applying external magnetic fields. For this pur¬ 
pose, we simulate a system with 8 x 8 x 512 spins includ¬ 
ing an easy-axis parallel to z-direction with dz = 10“^, 
a damping constant of a = 0.01 and apply an exter¬ 
nal field in the z-direction. Linear temperature gradients 
with variable lengths L excite net magnonic spin currents 
due to a non-equilibrium of the local magnonic density 
of states. In the hotter area more magnons exist than 
in the colder area and therefore more magnons propa¬ 
gate towards the colder region than the other way round, 
leading to a net magnon current from hot to cold. Since 
we aim to describe spin accumulation in samples with 
thin film geometry we will call the spatial extension of 
the gradient, L, from now on the thickness. 



FIG. 1 . Local magnon accumulation Am profile with the 
space position 2 in units of the lattice constant a for different 
thicknesses L. 


After an initial relaxation the system reaches a quasi 
static state where we calculate the local magnetization 
profile m{z) as an average over time and over the x-y- 
planes. The magnon accumulation can be defined as the 
difference of the local magnetization to its equilibrium 
value related to the local phonon temperature. 


Am{z) = m{z) - meq{Tp{z)). 


( 5 ) 


FIG.[l]shows exemplary the resulting magnon accumu¬ 
lation for an external field Bz = O.lJ/fis and a temper¬ 
ature gradient of VzT = 10~^J/{kBa) for various thick¬ 
nesses L. As shown in similar simulations by Kehlberger 
et al. the magnon accumulation has two extrema at the 
ends of the temperature gradient and their values in¬ 
crease with increasing thickness up to a characteristic 
length scale above which the magnon accumulation satu¬ 
rates. This characteristic length scale can be referred to 
as the mean magnon propagation length ■Cavg 01 ■ 

In the upper part of FIG. [21 the extreme value of the 
magnon accumulation at the cold end of the temperature 
gradient is shown versus applied magnetic field Bz for 
different thicknesses. The magnetic fields used in these 
simulations corresponds to magnetic fields of the order 
of 80 T leading to strong effects. In all cases one can see 
that the magnon accumulation decreases with increasing 
field, but the observed suppression of the magnon accu¬ 
mulation is thickness dependent. For thicker films, larger 
than the mean magnon propagation length, the field sup¬ 
pression of the magnon accumulation is stronger than for 
thinner films where the suppression effect shows only a 
weak thickness dependence. This can be seen in the lower 
part of FIG. | 2 ] where the normalized magnon accumula¬ 
tion Am{Bz)/Am{0) is shown for various thicknesses. 

The magnon accumulations at both, the hot and 
the cold end of the temperature gradient, are nearly 
proportional to the temperature difference between the 
magnonic and phononic subsystems 18 | . It was shown by 
Xiao et al. that this temperature difference scales linearly 
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magnetic field ^^Bz/ J 


FIG. 2. Upper part: Magnon accumulation at the cold end 
of the temperature gradient versus applied magnetic field Bz 
for different thicknesses with VT = lQ~^k^T/{Ja), a = 0.01, 
dz = 10~® J. Lower part: Corresponding normalized magnon 
accumulation Am{B)/Am{0). 


with the s pin Seebeck voltage, which is measured in ex¬ 
periments 17|. Hence, one can expect that the magnonic 
contribution to the SSE can be suppressed by magnetic 
fields. Interestingly, this behavior shows a thickness de¬ 
pendence with the influence of the field decreasing for 
thinner films, an effect that we will discuss later on. 

First, in order to verify our theoretical model, the high 
magnetic field dependence of the SSE is measured us¬ 
ing Pt/YIG hybrids. The samples used in the present 
study are (lll)-oriented YIG slabs (5x10x1 mm^) and 
YIG thin films grown on Gd 3 Ga 50 i 2 (GGG) substrate 
(5xl0x0.5mm^) by liquid phase epitaxy. The thickness 
of the YIG thin films is varied from 0.3 pm to 50 pm. 
5.5 nm-thick Pt layers are deposited on the YIG surfaces 
by dc-magnetron sputtering and further patterned into 
strips (length of 4 mm and width of 100 pm) by opti¬ 
cal lithography and argon ion beam milling. FIG. [Ua) 
shows a schematic diagram of the prepared sample and 
measurement setup. We adopt the longitudinal configu¬ 
ration to determine the SSE. The sample is sandwiched 
between the resistive heater and thermal sensor, further 
mounted onto a copper heat sink then put into the cryo¬ 
stat. Each functional layer is structured on a AI 2 O 3 sub¬ 
strate, which prevent an electrical short circuit the indi¬ 
vidual elements another, while ensuring maximal vertical 
temperature transport. Furthermore a thermal grease is 
used for the mounting of the elements, which ensures 
good thermal connection at the interfaces. A temper¬ 
ature gradient, VT, can be generated by the attached 
resistive heater and be varied by simply changing the 


heating currents. A spin current is thermally generated 
in the Pt layer along the direction of thermal gradient 
(^-axis), and further converted into an electric field due 
to the inverse spin Hall effect. The SSE can be detected 
electrically by measuring the voltage drop l^sE at the 
two ends of the Pt strip (along x-axis). The beneht 
of our measurement setup is that the Pt strips on the 
top and bottom surfaces of the samples can be utilized 
as an excellent resistance-temperature detector to deter¬ 
mine the temperature differences across the hybrids pre¬ 
cisely. In our room-temperature experiments we keep the 
cryostat at a fixed temperature of 300 K and a heating 
current in our resistive heater of /heat = 9 mA results 
in the temperature difference AT of roughly 4.2 ± 0.2 K 
for our measurements. An external magnetic field H is 
applied perpendicular to the Pt strip along y-axis. The 
maximum magnetic field is up to 9 T, which is orders 
of magnitude higher than the coercivity of YIG, ensur¬ 
ing that the magnetization of the YIG slab (or YIG thin 
films) is well aligned along the field direction. To exclude 
proximity effects from the Pt/YIG interface as the origin 
of the observed field suppression of the SSE, the magne¬ 
toresistance of the Pt layer is monitored, which shows no 
noticeable change beyond the saturation field of magne¬ 
tization of the YIG film. 

The VsSE is recorded as a function of magnetic fields. 
The SSE signal appears when a magnetic field is applied 
and flips its sign when the field reverses. We find the 
magnitude of the SSE signal suppressed by high magnetic 
fields at room temperature. The maximum value of Vsse 
is found at the smallest field interval of 0.2 T. The thick¬ 
ness dependent SSE coefficient ctsse = Vsse (0.2 T)/AT 
exhibits a similar trend to the one observed in our pre¬ 
vious work: the ctsse is enhanced with increasing film 
thickness and saturates above a characteristic length, 
demonstrating the bulk origin of magnonic spin current 
0. FIG. [31(b) shows examples of the normalized spin See¬ 
beck voltages Hsse(//)/Vsse(0.2 T) for selected Pt/YIG 
hybrids. Obviously, the suppression effect is reduced dra¬ 
matically with decreasing film thickness. The experimen¬ 
tal results are in agreement with the numerical results 
using our proposed theoretical model showing a field sup¬ 
pression up to 40% for the 1mm thick YIG slab under a 
magnetic field of 9T, while this value drops to only 0.1% 
for 0.3 |xm thick YIG films. 

To understand this effect at a fundamental level, we 
analyze it next based on an analytical linear spin wave 
theory. In particular we determine in this model the ori¬ 
gin of the held suppression of the magnon accumulation 
as well as the thickness dependence of this effect. The 
magnonic dispersion relation of the simulated system can 
be calculated by solving the linearized LEG equation, 

huj^i = 2dz + fisBz + 2J ^(1 - cosiqeae)) , (6) 

6 

where 6 denotes the spatial coordinates x,y,z. The 
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FIG. 4. Magnon frequency distribution of propagating 
magnons at the cold end of the temperature gradient and the 
frequency dependent magnon propagation length for different 
magnetic fields. 


The maximum propagation length is given by 


Cmax 

a 


j-.rz: 

\/2a V 


( 8 ) 


FIG. 3. (a) Schematic diagram of sample structure and 

measurement setup, (b) Normalized spin Seebeck voltage 
Vssb(4I)/Vsse(0.2 T) of Pt/YIG structures composed of YIG 
single crystal and selected YIG thin films with thickness rang¬ 
ing from 0.3 to 50 pm. 

dispersion relation consists of a magnon frequency gap 
?iWinin = + itgBz) and a second term depending on 

the wave vector q. The frequency gap increases with 
increasing magnetic fields and, hence, parts of the spec¬ 
trum are frozen out for higher magnetic fields. This effect 
can be seen in Fig. 21 By calculating a Fourier transfor¬ 
mation of the magnon accumulation Am in the time do¬ 
main, as used by Ritzmann et al. [isj ] , at the cold end of 
the temperature gradient, we can calculate the frequency 
distribution of the magnons from our simulations. The 
shown frequency distributions for temperature gradients 
with width L = 200a for different magnetic fields shows 
a clear shift of the frequency gap and therefore a sup¬ 
pression of parts of the frequency spectra. This effect 
can explain the suppression of the accumulation due to 
applied magnetic fields. 

The frequency dependent magnon propagation length 
^(w) is also shown in Fig. 2) It can be estimated using 
the lifetime r = and the group velocity of the 

magnons [3 leading to the expression: 



and can be much larger than the mean propagation 
length of the magnons. 

The thickness dependence of the field suppression can 
be explained by the frequency dependent magnon prop¬ 
agation length. Magnons with large propagation length 
will always reach the cold end of the temperature gradi¬ 
ent. When the thickness is reduced, the contribution of 
those magnons is reduced. Since only these magnons are 
suppressed by the frequency gap, the field suppression is 
also reduced for low thicknesses. Note, that the equa¬ 
tions above show that the same suppression effect would 
appear by modifying the anisotropy in the system, since 
the effect depends not only on the size of the magnetic 
field, but on the value of the frequency gap. 

Furthermore, not only parts of the spectrum are sup¬ 
pressed, but also the frequency dependent magnon prop¬ 
agation length ^(ui) is modified by applying higher mag¬ 
netic fields. When the propagation length of the magnons 
is reduced, the number of magnons propagating through 
the system is smaller, leading also to a reduction of the 
magnon accumulation, showing that multiple effects con¬ 
tribute to this generic property of thermal magnon prop¬ 
agation in applied magnetic fields. 

In conclusion we have shown numerical simulations and 
experiments on the magnetic field suppression of the lon¬ 
gitudinal SSE. Applying large magnetic fields can sup¬ 
press the excited magnon spin currents by suppressing 
parts of the frequency spectrum. We find qualitative 
agreement between the performed simulation and exper¬ 
iments of the field dependence of the longitudinal SSE in 
YIG at room temperature for various thicknesses. The 
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suppression is strongly dependent on the thickness of the 
sample and is practically vanishing for thin films since low 
frequency magnons play a minor role in thin films as com¬ 
pared to bulk materials. Furthermore, not only parts of 
the spectra are suppressed but the magnon propagation 
length given by the Gilbert damping is modified lead¬ 
ing to a lower propagation length in systems with higher 
magnetic fields. Our analysis opens a new avenue to de¬ 
termining which magnons contribute to the spin current 
as a careful study of the field dependence would allow 
then to measure the intensities of the involved frequency 
as proposed by Boona et al. 2l|. Furthermore, the held 
suppression open new opportunities for the control of the 
magnon flow and resulting SSE with external helds. Dur¬ 
ing the preparation of this manuscript we became aware 
of two related experimental works in arxiv 25, 2^ sup¬ 
porting our hndings. 
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